Persistence of neutral polymorphisms in Lake Victoria cichlid fish

Phylogenetic trees for groups of closely related species often have different topologies, depending on the genes used. One explanation for the discordant topologies is the persistence of polymorphisms through the speciation phase, followed by differential fixation of alleles in the resulting species. The existence of transspecies polymorphisms has been documented for alleles maintained by balancing selection but not for neutral alleles. In the present study, transspecific persistence of neutral polymorphisms was tested in the endemic haplochromine species flock of Lake Victoria cichlid fish. Putative noncoding region polymorphisms were identified at four randomly selected nuclear loci and tested on a collection of 12 Lake Victoria species and their putative riverine ancestors. At all loci, the same polymorphism was found to be present in nearly all the tested species, both lacustrine and riverine. Different polymorphisms at these loci were found in cichlids of other East African lakes (Malawi and Tanganyika). The Lake Victoria polymorphisms must have therefore arisen after the flocks now inhabiting the three great lakes diverged from one another, but before the riverine ancestors of the Lake Victoria flock colonized the Lake. Calculations based on the mtDNA clock suggest that the polymorphisms have persisted for about 1.4 million years. To maintain neutral polymorphisms for such a long time, the population size must have remained large throughout the entire period.

Methylenetetrahydrofolate reductase (MTHFR) polymorphisms and risk of molecularly defined subtypes of childhood acute leukemia

Low folate intake as well as alterations in folate metabolism as a result of polymorphisms in the enzyme methylenetetrahydrofolate reductase (MTHFR) have been associated with an increased incidence of neural tube defects, vascular disease, and some cancers. Polymorphic variants of MTHFR lead to enhanced thymidine pools and better quality DNA synthesis that could afford some protection from the development of leukemias, particularly those with translocations. We now report associations of MTHFR polymorphisms in three subgroups of pediatric leukemias: infant lymphoblastic or myeloblastic leukemias with MLL rearrangements and childhood lymphoblastic leukemias with either TEL-AML1 fusions or hyperdiploid karyotypes. Pediatric leukemia patients (n = 253 total) and healthy newborn controls (n = 200) were genotyped for MTHFR polymorphisms at nucleotides 677 (C→T) and 1,298 (A→C). A significant association for carriers of C677T was demonstrated for leukemias with MLL translocations (MLL+, n = 37) when compared with controls [adjusted odd ratios (OR) = 0.36 with a 95% confidence interval (CI) of 0.15–0.85; P = 0.017]. This protective effect was not evident for A1298C alleles (OR = 1.14). In contrast, associations for A1298C homozygotes (CC; OR = 0.26 with a 95% CI of 0.07–0.81) and C677T homozygotes (TT; OR = 0.49 with a 95% CI of 0.20–1.17) were observed for hyperdiploid leukemias (n = 138). No significant associations were evident for either polymorphism with TEL-AML1+ leukemias (n = 78). These differences in allelic associations may point to discrete attributes of the two alleles in their ability to alter folate and one-carbon metabolite pools and impact after DNA synthesis and methylation pathways, but should be viewed cautiously pending larger follow-up studies. The data provide evidence that molecularly defined subgroups of pediatric leukemias have different etiologies and also suggest a role of folate in the development of childhood leukemia.

Asian genotypes of JC virus in Native Americans and in a Pacific Island population: Markers of viral evolution and human migration

The human polyomavirus JC (JCV) causes the central nervous system demyelinating disease progressive multifocal leukoencephalopathy. Previously, we showed that 40% of Caucasians in the United States excrete JCV in the urine as detected by PCR. We have now studied 68 Navaho from New Mexico, 25 Flathead from Montana, and 29 Chamorro from Guam. By using PCR amplification of a fragment of the VP1 gene, JCV DNA was detected in the urine of 45 (66%) Navaho, 14 (56%) Flathead, and 20 (69%) Chamorro. Genotyping of viral DNAs in these cohorts by cycle sequencing showed predominantly type 2 (Asian), rather than type 1 (European). Type 1 is the major type in the United States and Hungary. Type 2 can be further subdivided into 2A, 2B, and 2C. Type 2A is found in China and Japan. Type 2B is a subtype related to the East Asian type, and is now found in Europe and the United States. The large majority (56–89%) of strains excreted by Native Americans and Pacific Islanders were the type 2A subtype, consistent with the origin of these strains in Asia. These findings indicate that JCV infection of Native Americans predates contact with Europeans, and likely predates migration of Amerind ancestors across the Bering land bridge around 12,000–30,000 years ago. If JCV had already differentiated into stable modern genotypes and subtypes prior to first settlement, the origin of JCV in humans may date from 50,000 to 100,000 years ago or more. We conclude that JCV may have coevolved with the human species, and that it provides a convenient marker for human migrations in both prehistoric and modern times.

Direct genetic analysis by matrix-assisted laser desorption/ionization mass spectrometry

An approach to analyzing single-nucleotide polymorphisms (SNPs) found in the human genome has been developed that couples a recently developed invasive cleavage assay for nucleic acids with detection by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). The invasive cleavage assay is a signal amplification method that enables the analysis of SNPs by MALDI-TOF MS directly from human genomic DNA without the need for initial target amplification by PCR. The results presented here show the successful genotyping by this approach of twelve SNPs located randomly throughout the human genome. Conventional Sanger sequencing of these SNP positions confirmed the accuracy of the MALDI-TOF MS analysis results. The ability to unambiguously detect both homozygous and heterozygous genotypes is clearly demonstrated. The elimination of the need for target amplification by PCR, combined with the inherently rapid and accurate nature of detection by MALDI-TOF MS, gives this approach unique and significant advantages in the high-throughput genotyping of large numbers of SNPs, useful for locating, identifying, and characterizing the function of specific genes.

The major murine systemic lupus erythematosus susceptibility locus, Sle1, is a cluster of functionally related genes

The major murine systemic lupus erythematosus (SLE) susceptibility locus Sle1 is syntenic to a chromosomal region linked with SLE susceptibility in multiple human studies. Congenic analyses have shown that Sle1 breaks tolerance to chromatin, a necessary step for full disease induction that can be suppressed by specific modifier loci. In the present study, our fine mapping analysis of the location of Sle1 has determined that three loci within this congenic interval, termed Sle1a, Sle1b, and Sle1c, can independently cause a loss of tolerance to chromatin. Each displays a distinctive profile of serological and cellular characteristics, with T and B cell functions being more affected by Sle1a and Sle1b, respectively. The epistatic interactions of Sle1 with other susceptibility loci to cause severe nephritis cannot be accounted, however, by these three loci alone, suggesting the existence of an additional locus, termed Sle1d. These findings indicate that the potent autoimmune phenotype caused by the Sle1 genomic interval reflects the combined impact of four, separate, susceptibility genes. This level of genetic complexity, combined with similar findings in other systems, supports the possibility that many complex trait loci reflect the impact of polymorphisms in linked clusters of genes with related functions.

Characterization of repetitive DNA in the Mycoplasma genitalium genome: possible role in the generation of antigenic variation.

We have characterized a family of repetitive DNA elements with homology to the MgPa cellular adhesion operon of Mycoplasma genitalium, a bacterium that has the smallest known genome of any free-living organism. One element, 2272 bp in length and flanked by DNA with no homology to MgPa, was completely sequenced. At least four others were partially sequenced. The complete element is a composite of six regions. Five of these regions show sequence similarity with nonadjacent segments of genes of the MgPa operon. The sixth region, located near the center of the element, is an A+T-rich sequence that has only been found in this repeat family. Open reading frames are present within the five individual regions showing sequence homology to MgPa and the adjacent open reading frame 3 (ORF3) gene. However, termination codons are found between adjacent regions of homology to the MgPa operon and in the A+T-rich sequence. Thus, these repetitive elements do not appear to be directly expressible protein coding sequences. The sequence of one region from five different repetitive elements was compared with the homologous region of the MgPa gene from the type strain G37 and four newly isolated M. genitalium strains. Recombination between repetitive elements of strain G37 and the MgPa operon can explain the majority of polymorphisms within our partial sequences of the MgPa genes of the new isolates. Therefore, we propose that the repetitive elements of M. genitalium provide a reservoir of sequence that contributes to antigenic variation in proteins of the MgPa cellular adhesion operon.

HLA-Cw allele analysis by PCR-restriction fragment length polymorphism: study of known and additional alleles.

We describe a technique for HLA-Cw genotyping by digestion of PCR-amplified genes with restriction endonucleases. Locus-specific primers selectively amplified HLA-Cw sequences from exon 2 in a single PCR that avoided coamplification of other classical and nonclassical class I genes. Amplified DNAs were digested with selected enzymes. Sixty-three homozygous cell lines from International Histocompatibility Workshop X and 113 unrelated individual cells were genotypes for HLA-Cw and compared with serology. The present protocol can distinguish 23 alleles corresponding to the known HLA-Cw sequences. Genotyping of serologically undetectable alleles (HLA-Cw Blank) and of heterozygous cells was made possible by using this method. Six additional HLA-Cw alleles were identified by unusual restriction patterns and confirmed by sequencing; this observation suggests the presence of another family of allele-sharing clusters in the HLA-B locus. This PCR-restriction endonuclease method provides a simple and convenient approach for HLA-Cw DNA typing, allowing the definition of serologically undetectable alleles, and will contribute to the evaluation of the biological role of the HLA-C locus.

Mouse mammary tumor virus/v-Ha-ras transgene-induced mammary tumors exhibit strain-specific allelic loss on mouse chromosome 4

Hybrid mice carrying oncogenic transgenes afford powerful systems for investigating loss of heterozygosity (LOH) in tumors. Here, we apply this approach to a neoplasm of key importance in human medicine: mammary carcinoma. We performed a whole genome search for LOH using the mouse mammary tumor virus/v-Ha-ras mammary carcinoma model in female (FVB/N × Mus musculus castaneus)F1 mice. Mammary tumors developed as expected, as well as a few tumors of a second type (uterine leiomyosarcoma) not previously associated with this transgene. Genotyping of 94 anatomically independent tumors revealed high-frequency LOH (≈38%) for markers on chromosome 4. A marked allelic bias was observed, with M. musculus castaneus alleles almost exclusively being lost. No evidence of genomic imprinting effects was noted. These data point to the presence of a tumor suppressor gene(s) on mouse chromosome 4 involved in mammary carcinogenesis induced by mutant H-ras expression, and for which a significant functional difference may exist between the M. musculus castaneus and FVB/N alleles. Provisional subchromosomal localization of this gene, designated Loh-3, can be made to a distal segment having syntenic correspondence to human chromosome 1p; LOH in this latter region is observed in several human malignancies, including breast cancers. Evidence was also obtained for a possible second locus associated with LOH with less marked allele bias on proximal chromosome 4.

A novel modifier gene for plasma von Willebrand factor level maps to distal mouse chromosome 11

Type 1 von Willebrand disease (VWD), characterized by reduced levels of plasma von Willebrand factor (VWF), is the most common inherited bleeding disorder in humans. Penetrance of VWD is incomplete, and expression of the bleeding phenotype is highly variable. In addition, plasma VWF levels vary widely among normal individuals. To identify genes that influence VWF level, we analyzed a genetic cross between RIIIS/J and CASA/Rk, two strains of mice that exhibit a 20-fold difference in plasma VWF level. DNA samples from F2 progeny demonstrating either extremely high or extremely low plasma VWF levels were pooled and genotyped for 41 markers spanning the autosomal genome. A novel locus accounting for 63% of the total variance in VWF level was mapped to distal mouse chromosome 11, which is distinct from the murine Vwf locus on chromosome 6. We designated this locus Mvwf for “modifier of VWF.” Additional genotyping of as many as 2407 meioses established a high resolution genetic map with gene order Cola1-Itg3a-Ngfr-Mvwf/Gip-Hoxb9-Hoxb1-Cbx·rs2-Cox5a-Gfap. The Mvwf candidate interval between Ngfr and Hoxb9 is ≈0.5 centimorgan (cM). These results demonstrate that a single dominant gene accounts for the low VWF phenotype of RIIIS/J mice in crosses with several other strains. The pattern of inheritance suggests a gain-of-function mutation in a unique component of VWF biosynthesis or processing. Characterization of the human homologue for Mvwf may have relevance for a subset of type 1 VWD cases and may define an important genetic factor modifying penetrance and expression of mutations at the VWF locus.

Direct sequencing of bacterial and P1 artificial chromosome-nested deletions for identifying position-specific single-nucleotide polymorphisms

A loxP-transposon retrofitting strategy for generating large nested deletions from one end of the insert DNA in bacterial artificial chromosomes and P1 artificial chromosomes was described recently [Chatterjee, P. K. & Coren, J. S. (1997) Nucleic Acids Res. 25, 2205–2212]. In this report, we combine this procedure with direct sequencing of nested-deletion templates by using primers located in the transposon end to illustrate its value for position-specific single-nucleotide polymorphism (SNP) discovery from chosen regions of large insert clones. A simple ampicillin sensitivity screen was developed to facilitate identification and recovery of deletion clones free of transduced transposon plasmid. This directed approach requires minimal DNA sequencing, and no in vitro subclone library generation; positionally oriented SNPs are a consequence of the method. The procedure is used to discover new SNPs as well as physically map those identified from random subcloned libraries or sequence databases. The deletion templates, positioned SNPs, and markers are also used to orient large insert clones into a contig. The deletion clone can serve as a ready resource for future functional genomic studies because each carries a mammalian cell-specific antibiotic resistance gene from the transposon. Furthermore, the technique should be especially applicable to the analysis of genomes for which a full genome sequence or radiation hybrid cell lines are unavailable.

Scrapie susceptibility-linked polymorphisms modulate the in vitro conversion of sheep prion protein to protease-resistant forms

Prion diseases are natural transmissible neurodegenerative disorders in humans and animals. They are characterized by the accumulation of a protease-resistant scrapie-associated prion protein (PrPSc) of the host-encoded cellular prion protein (PrPC) mainly in the central nervous system. Polymorphisms in the PrP gene are linked to differences in susceptibility for prion diseases. The mechanisms underlying these effects are still unknown. Here we describe studies of the influence of sheep PrP polymorphisms on the conversion of PrPC into protease-resistant forms. In a cell-free system, sheep PrPSc induced the conversion of sheep PrPC into protease-resistant PrP (PrP-res) similar or identical to PrPSc. Polymorphisms present in either PrPC or PrPSc had dramatic effects on the cell-free conversion efficiencies. The PrP variant associated with a high susceptibility to scrapie and short survival times of scrapie-affected sheep was efficiently converted into PrP-res. The wild-type PrP variant associated with a neutral effect on susceptibility and intermediate survival times was converted with intermediate efficiency. The PrP variant associated with scrapie resistance and long survival times was poorly converted. Thus the in vitro conversion characteristics of the sheep PrP variants reflect their linkage with scrapie susceptibility and survival times of scrapie-affected sheep. The modulating effect of the polymorphisms in PrPC and PrPSc on the cell-free conversion characteristics suggests that, besides the species barrier, polymorphism barriers play a significant role in the transmissibility of prion diseases.

Genetic epidemiology of single-nucleotide polymorphisms

On the causal hypothesis, most genetic determinants of disease are single-nucleotide polymorphisms (SNPs) that are likely to be selected as markers for positional cloning. On the proximity hypothesis, most disease determinants will not be included among markers but may be detected through linkage disequilibrium with other SNPs. In that event, allelic association among SNPs is an essential factor in positional cloning. Recent simulation based on monotonic population expansion suggests that useful association does not usually extend beyond 3 kb. This is contradicted by significant disequilibrium at much greater distances, with corresponding reduction in the number of SNPs required for a cost-effective genome scan. A plausible explanation is that cyclical expansions follow population bottlenecks that establish new disequilibria. Data on more than 1,000 locus pairs indicate that most disequilibria trace to the Neolithic, with no apparent difference between haplotypes that are random or selected through a major disease gene. Short duration may be characteristic of alleles contributing to disease susceptibility and haplotypes characteristic of particular ethnic groups. Alleles that are highly polymorphic in all ethnic groups may be older, neutral, or advantageous, in weak disequilibrium with nearby markers, and therefore less useful for positional cloning of disease genes. Significant disequilibrium at large distance makes the number of suitably chosen SNPs required for genome screening as small as 30,000, or 1 per 100 kb, with greater density (including less common SNPs) reserved for candidate regions.

Sensitive detection of DNA polymorphisms by the serial invasive signal amplification reaction

The invasive signal amplification reaction has been previously developed for quantitative detection of nucleic acids and discrimination of single-nucleotide polymorphisms. Here we describe a method that couples two invasive reactions into a serial isothermal homogeneous assay using fluorescence resonance energy transfer detection. The serial version of the assay generates more than 107 reporter molecules for each molecule of target DNA in a 4-h reaction; this sensitivity, coupled with the exquisite specificity of the reaction, is sufficient for direct detection of less than 1,000 target molecules with no prior target amplification. Here we present a kinetic analysis of the parameters affecting signal and background generation in the serial invasive signal amplification reaction and describe a simple kinetic model of the assay. We demonstrate the ability of the assay to detect as few as 600 copies of the methylene tetrahydrofolate reductase gene in samples of human genomic DNA. We also demonstrate the ability of the assay to discriminate single base differences in this gene by using 20 ng of human genomic DNA.